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■ We complete the analysis of twist-two generalized parton distributions of the nucleon in 

one-loop order of heavy-baryon chiral perturbation theory. Extending our previous study 
Qh! of the chiral-even isosinglet sector, we give results for chiral-even isotriplet distributions 
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^ ■ generalized parton distributions of the pion. 
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1 Introduction 



Generalized parton distributions (GPDs) provide a unified parameterization of many different aspects 
of liadron physics [U EJ [3l Hj • Understanding GPDs in detail is therefore tantamount to understanding 
in large parts the internal structure of hadrons. This motivates extensive experimental programs 
as well as theoretical work. Details can be found in the reviews O El [7], which emphasize the 
different types of physics encoded in these quantities. More recently it has been shown that interesting 
information about the distribution of transversely polarized quarks in a hadron is contained in GPDs 
associated with chiral-odd quark operators [51 [9], for which there have been relatively few studies so 
far. 

The extraction of GPDs from experiment is a highly non-trivial task, since in observables the 
distributions appear only within convolutions. These are relatively simple at leading order in the 
strong coupling but become increasingly complex at higher orders, see e.g. [lOj . In practice one 
therefore has to use parameterizations of GPDs which are on one hand sufficiently flexible to catch 
the physics and on the other hand contain only few parameters. In this context, the calculation of 
moments of GPDs in lattice QCD is expected to become highly important in the future. 

The lattice evaluation of these moments, parameterized by the form factors of local matrix ele- 
ments, is very similar to the case of the usual electromagnetic form factors [llj. The main limitation 
at present is that lattice calculations with dynamical quarks can only be done for unphysically heavy 
quarks and thus pions. The mass of the pion affects however the spatial extent of the nucleon and 
hence its form factors. Therefore, their extrapolation to the physical limit can be fairly non-trivial, 
and simple linear extrapolations with respect to ni^ or could be quite inadequate. Progress in 
this respect requires an analysis within chiral perturbation theory (ChPT). We have presented such 
an analysis for the pion GPDs in [12] and for nucleon GPDs in the chiral-even isosinglet sector in [13]. 
In the present paper we extend this work to the chiral-even isotriplet sector and the chiral-odd sector, 
giving complete corrections at one-loop accuracy. Calculations of a similar scope have recently been 
reported in [T3], and we will compare our results in detail. There already exists a number of lattice 
results for moments of GPDs, see |151 ITB] and references therein. We do not include any ChPT fits 
to these in the present paper, but leave them to future lattice studies. 

Our paper is organized as follows. In Sections O [3] and [6] we collect details about GPD parameter- 
izations, the operator product expansion, and heavy-baryon ChPT that are needed in our analysis. 
We proceed in each case by constructing the operators within ChPT that match the relevant twist-two 
operators in QCD, and by identifying the loop corrections which contribute to a given form factor 
at relative order 0{q^) in the chiral expansion (Sections [H 15.11 and 171). Results of the corresponding 
calculations are given for the vector form factors in Section 15.21 for the axial form factors in Section 
15.31 and for the chiral-odd form factors in Section [8l In Section [9] all results are collected and rewritten 
in terms of the usual parameterization of GPDs. We summarize our main findings in Section [TOl 

2 Chiral-even generalized parton distributions 

To begin with let us recall the definitions of generalized parton distributions associated with chiral- 
even quark operators. For the distributions with definite isospin / in a nucleon one can write 
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where a is a light-like auxiliary vector, M is the nucleon mass, and we use the standard kinematical 
variables P = ^{p + p'), A = p' — p, t = and 2^ = — (Aa)/(Pa). Wilson lines must be inserted 
between the quark fields if one is not working in the light-cone gauge (aA) = 0. We combine the 
two-dimensional unit matrix r*^ and the triplet of Pauli matrices r in a four- vector = (r*^,-?), 
with the matrices acting on the isodoublet of quark fields q or of nucleon states N. The isosinglet 
distributions correspond to A = and the isotriplet ones to ^ = 1, 2, 3. In terms of individual quark 
flavors in the proton one has H^^^ = -\- H'^ and H^^^ = — H'^, with analogous relations for 
the other distributions. 

The Mellin moments of the GPDs in ([T|) are related to the matrix elements of the chiral-even local 
twist-two operators 
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with = ^(-D^ — -D^). Here T denotes the subtraction of trace terms in the indicated Lorentz 
indices and S denotes symmetrization, normalized as S^^^j t^^^'^ = ^{t^^^^ +t^^'^'^). Both operations 
are conveniently implemented by contraction with the auxiliary vector a, 

The local matrix elements can be parameterized as 
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and the moments of the GPDs are given by 



dx H{x, t) = Y,{2e,f An,k{t) + mod(n + 1, 2) (2C)"C7„(t) 

even 

»1 n-l 

/ dx x"-i E{x, t) = 5^(20' - mod(n + 1, 2) (20"C„(t) 

J —I I — n 



fc=0 
even 

n-l 



^ dx x''~^ Hix, t) = ^(20' A„,fc(t) , 



k=0 
even 

n-l 



/ dxx''-^E{x,C,t) = y2i2C)''Bn,k{t), 

J-l tTn 



(5) 



fc=0 
even 



3 



where here and in the following we omit the isospin label / when it is not required. The restriction 
to even A: in (jj]) and ([5]) is a consequence of time reversal invar iance. 

To calculate the chiral corrections to the nucleon form factors in heavy-baryon chiral perturbation 
theory we work in the Breit frame, where P = 0. The incoming and outgoing nucleons then have 
opposite spatial momenta p' = —p = A/2 and equal energies, Pq = Pq = M7 with 



7 = Vl- A2/4M2 



(6) 



In terms of the velocity vector v, given hy v = (1, 0, 0, 0) in the Breit frame, the incoming and outgoing 
nucleon momenta are given hy p = M'jv — A/2 and p' = M'yv + A/2. Note that (vA) = (vS) = 0. 
Dirac bilinears can be expressed in terms of the velocity and the spin operator = ^icTf^^'jc^v'^ . 
Introducing the spinors 



Uv{p) = A/" 



-1 



^u{p), u,(p')=M-'^u{p') 



(7) 



with J\f = y/ {1 + 7)/2, the matrix elements in Q can be rewritten as [13] 

n-l 
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where due to time reversal invariance the terms with E^.k+i are only nonzero for odd A;, whereas 
those with M^^k^ E^^k and M„ ^ are only nonzero for even k. The relation between the form factors 
in dH) and those in dHl) is 



A2 
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which is readily inverted to 
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3 Heavy-baryon ChPT 

To set our notation, let us briefly review the main ingredients of chiral perturbation theory for heavy 
baryons, which is an effective theory for the limit q, <C M, where g is a generic momentum. 
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To describe pions we use the nonlinear representation Uix) = [n(x)]^ = exp [i7r"(x) t"/-F] , where 
F ~ 92 MeV is the pion decay constant in the chiral hmit|l| The exphcit breaking of chiral symmetry 
by the quark masses is implemented by the field x(^)- We assume the isospin limit, where one can 
replace x{x) ~^ w?t^ with the bare pion mass m. We will not use external vector or axial vector 
fields here. The nucleon is described by the heavy-baryon field Ny{x) = ^(1 + f) e**^"*-^^^ A^(x), where 
Mq is the bare nucleon mass and v the velocity vector. The Fourier transform of N^{x) depends on 
the residual nucleon momentum, given by the original nucleon momentum minus Mqv. Important 
derived quantities are the axial vector field 

Uf, = i{u^df,u - udf^u^) = S^vrV" + 0(7r3) , (11) 



the connection 



and 



= ^ (u^d^u + ud^u^^ = e-"^ d^n'r^ + 0{7r^) , (12) 



X± = u^Xu^ ± ux^u . (13) 



Under global chiral transformations, described by unitary matrices Vl and Vr, the different fields 
transform as 



u 



r^^HT^H^ + Hd^H\ (14) 

and and x± transform homogeneously as 

u^^Hu^H^, X±^Hx±H^. (15) 

The unitary matrix H depends on Vl, and on U{x) and therefore has an x dependence. With 
the connection one can construct the covariant derivative V^. It acts as V^X = d^X + T^X on 
quantities like A''^,, which transform with a factor H on their left, and as [V^,y] = dfj^Y + [r^,y] 
on quantities like u^, which transform with H on the left and with on the right. Corresponding 
derivatives acting to the left are ZV = Zd — ZT^ and V] = Yd — [y, F^], where Z transforms 
with a factor H'^ on its right. 

The effective Lagrangian for the theory contains a pure pion piece and a piece describing the 
nucleon and its interaction with pions, C^s = ^tt + C--kN- Expanding in powers of q one has 

with [m [S] 

42) = :^Tr(n^n^ + X+), 

= i| {Trx+Y + ^ {2Trx+ Triu^u^ +2TTix'_) - (Tr X-)'} + • • • , (17) 



'^Our convention is that uppercase indices of r as in ([1} run from to 3, whereas lowercase ones run from 1 to 3. 
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and 119] 



6^^ -N 



'^{^i{vV)+go(,Su)^N^, 



+ (c2 - ^) (^^n)' + C3 n,n^ + (c4 + ^) 51 n^n.j TV^ , (18) 

where go is the nucleon axial- vector coupling in the chiral hmit and the li and q are further low-energy 
constants. The terms not displayed in C^t^^ couple to at least four pion fields and will not be needed 
in our calculations. 

For calculating nucleon matrix elements in the Breit frame we need the residual momenta of the 
incoming and outgoing nucleon, 

r = p - Mqv = wv - A/2 , r' = p - Mqv = wv + A/2 (19) 

with 

yj = M{j -1) +5M = -— -Acim'^ + 0{q^) , (20) 

oM 

where 5M = M — Mq is the nucleon mass shift. Using the spinors ([7]) one obtains a matrix element 
as [20] 

{p'\0\p)=M^ZNUv{p)Go{r',r)u,{p), (21) 

where Go{r' ,r) is the truncated Green function for external heavy-baryon fields A^^,, and the 
operator O in the effective theory. Zn is the heavy-baryon field renormalization constant, 

= 1 - 2IW - t - + ■ '''' 

where d2^{ii) is a low-energy constant in the Lagrangian C^j^ given in [21]. 

4 Chiral even isotriplet operators 
4.1 Construction of effective operators 

To find the operators in the effective theory which match the quark-gluon operators ([2|) in QCD we 
generalize the construction of [13] to the isotriplet sector. The relevant effective operators contain a 
part which involves only pion fields (and couples to the nucleon via interactions from C-j^n) and 
a part OttA^ that is bilinear in the nucleon field. We thus have 

0^{a) - OlM + O^Nia) , d^{a) - O^^M + On,.iv(a) , (23) 

where for the pure pion operators 0^^{a) and 0^^{a) we will use the form given in [12j. The 

pion-nucleon operators 0^^^{a) and 0^^^{a) are conveniently constructed by first matching the 
operators 

(0^(a))^, = i {laDT-^ q, , (0^(a))„. = q^ </ (iaB)""! q, , (24) 
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where i and j are isospin indices. They involve quarks of definite chirahty and transform as 

0^{a) ^ V^O^{a)V^ 0^{a) ^ V^O^{a)vl (25) 

unter chiral rotations. Parity transforms 0^{a) and Oi^{a) into each other. The corresponding 
effective operators that are bihnear in the nucleon field can be written in the form 

{Q^{a)).^ = (iVLOint)^. {u02N,)^ , (Q^(a)),^. = {N,0[u)^ {u^O'.N,). , (26) 

where Oi , O2 transform like under chiral rotations and O'l , O2 are related to them by parity. The 
vector and axial vector operators are then readily obtained as 

Ol^^{a) = Tvr^{Q^{a) + Q^(a)} Oi^Nia) = Trr^{Q^(a) - Q^(a)} (27) 

and will involve the combinations 

rf^ = nV^n ± -ur^n^ , (28) 

where the subscript e indicates that they occur in chiral even operators. In the isosinglet case one 
has simply Tg^. = 2t^ and = 0, whereas the isotriplet combinations 

= 2r'^ + t^Htt^t^ - ^V) + 0{ti^) <_ = -|; e^^^nK" + 0(^3) (29) 

F F 

involve an even or odd number of pion fields, respectively. The operators Oi, O2 can be constructed 
from the fields and x±i from the covariant derivatives introduced in Section [3l One can 
rearrange the covariant derivatives in Q^{a) and Qn{a) to act either as total derivatives 9^ on the 
product of all fields or in the antisymmetric form V^j = ^(V^ — V^), where = d ^ + T ^ acts on 
the product of all fields to the right and = 9^ — T^ on the product of all fields to the left. The 
operators Q^(a) and Q^(a) are tensors having n indices contracted with the auxiliary vector a. Other 
than d^, and these tensors can contain the vectors and and the totally antisymmetric 
tensor. The number of spin vectors can be changed using the identities 

1 i 
{S\, S^} = -{vxVf, - gx^) , [S'a, S^] = iex^^up Sx = -- €x^,up v^iS", S^] (30) 

where our convention for the totally antisymmetric tensor is £0123 ~ ^- -^o^ operators under 
discussion we chose a basis where appears at most linearly, or quadratically as the commutator 
[^A, Sp]. For counting powers of q one associates chiral dimension 1 to d^, V^, and chiral dimension 
2 to x±- 

We now make all factors of (av) explicit and write 

n 

= M^-^-' {avT-^ Oi,{a) , (31) 

fc=0 

where 0^^,{a) is free of factors [av). For contracting the k vectors in 0^j^{a) one can use 
only once, so that this operator contains at least k — 1 vectors d^, or u^. Thus we can further 
decompose 

00 

On,fe(«)= (32) 
i=-l 

where 0^j^-{a) has chiral dimension k + i. For 0,^^^(a) one has a decomposition in full analogy to 
(i3T]) and ft. 
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Table 1: Overview of contributions to the chiral even form factors. The restriction in the second 
column is due to time reversal invariance. A^a is the number of factors (aA) and {SA) in the 
decomposition ([8]). The indices of the operators must satisfy I > k and i > 0, and the corresponding 
graphs contribute to the form factor at order 0{q'^) with d > — A^a and from ()33p . 



form factor 


k 


A^A 


operators 


En,k+1 


odd 


k + 1 


On,l+l,i-l 


Mn,k 


even 


k + 1 


On,l+l,i-l 




even 


k 


On,l+l,i-l 


Mn,k 


even 


k + 2 


On,l+l,i-l 



4.2 Power counting for tree and loop graphs 

As shown in |13j the chiral dimension of a graph with two external nucleon legs and insertion of the 
operator 0^^^^^{a) or O^^^ia) is 

Dk^, = 2L + k + i + Y^ (dimy^(j) - 2) + ^ (dimy^7v(j) - l) , (33) 

where L is the number of loops (with L = for tree graphs). Kr(j) and T4-Ar(j) respectively denote 
the jth vertex from C-,^ and Ct^n in the graph, A^^^- and Nj^n are the corresponding total numbers of 
vertices, and and Ij^j are the numbers of pion and nucleon propagators. Corrections to the nucleon 
propagator from higher orders of C^^n are counted as a nucleon-nucleon vertex and are accompanied 
by two (leading-order) nucleon propagators on either side. Notice that + r'^ = 2wv^ is of order 
0{q^) and thus one order higher than the generic power associated with a residual nucleon momentum. 
A graph with chiral dimension j can thus generate contributions to a nucleon matrix element of 
order 0{q'^) with d > D^ i. Since 0^f^^{a) is accompanied by a factor (af)""*' in ([3T]) it can only 
contribute to form factors with at least n — k powers of {av) in the decomposition ([8]) of the nucleon 
matrix element. Taking into account the number A'^a of factors (aA) and (SA) in that decomposition, 
one can establish the order in the chiral expansion to which a given operator can contribute to a form 
factor. The result is given in Table [H 

Throughout this paper we refer to orders 0{q'^) in the chiral expansion of a given form factor 
rather than the expansion of the corresponding matrix element. This is most convenient for the 
problem at hand, since the chiral order of matrix elements increases with the order n of the operator, 
whereas the chiral order of the form factors has as a natural point of reference the order 0{q^) from 
tree-level insertions of operators with the lowest chiral dimension at given n. 

The contributions of the operators C'^fc j(a) and 0:^i^-{a) at tree level are readily evaluated. The 
tree level graphs do not contain pions, so that one can replace 

^0, ^ iA^, V'^ ^ -iwv^, 

t^+^2t^, rt^O, x+^2mV, X- ^ 0. (34) 

Operators with do not contribute to the form factors at leading order since w is of order 0{q^). 
The different types of higher-order contributions to the form factors are discussed in Section 3.2 of 
|13j . In the results we give for the form factors, we lump them all into coefficients describing the 
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a be 

Figure 1: One-loop graphs with the insertion of a pion-nucleon operator On^nNio) or On,-KN{o)-, which 
is denoted by a black blob. Not shown is the analog of graph c with residual momentum / + wv + A/2 
of the intermediate nucleon line. 



and t corrections from tree graphs, except for the terms proportional to Qq in the expression (|22p of 
the wave function renormalization constant Z^v, which we combine with the terms due to loop graphs. 

The one-loop graphs with pion-nucleon operator insertions are shown in Fig. [TJ The construction 
of operators detailed in Section [4.11 allows one to easily track the origin of factors which arise from 
a graph and must match the factors in the form factor decomposition ([5]). For this we use that the 
denominators of the pion and nucleon propagators are — + and {Iv + w + respectively, 
so that the loop integration turns tensors l^^ ■ ■ ■ Ifj.^ into tensors constructed from and gxfj_. We find 
that with the leading-order (LO) interactions from C^jf and the next-to-leading (NLO) interactions 
from /^2v ^ factor A^ which is not contracted to A^ (and hence can be contracted with or S^^) 
can only originate from [13] 

1. a total derivative in the operator insertion, 

2. a term {lv){SA) due to an NLO pion-nucleon vertex, 

3. or a term (/A) due to an NLO nucleon propagator correction. 

We further find that two factors of A^ which are not contracted to A^ can originate as (/A) (AS) 
from the NNLO pion-nucleon vertex generated by the term 

- ^ iV; {(V5)(nV) + (Vn)(5V)} N, (35) 

(3) 

of the Lagrangian given in li21j. 

In the next section we will see that i-i(^) ^ni+i i-ii'^) with i = 0, 1, 2 have at most 

I + i total derivatives d^. A one-loop graph with insertion of such an operator and pion-nucleon 
interactions up to NNLO must therefore satisfy 

Z + i + ^ (dim Ktv (i) - 1) > A^A (36) 
i=i 

in order to produce the number of factors A^ required to contribute to the form factors in Table [H 
For i > 2 and for pion-nucleon interactions higher than NNLO this inequality is trivially fulfilled. 
With the power counting established in the table, one then finds that the one-loop contributions from 
pion-nucleon operators for all form factors start at order 0{(f'). 
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a be 

Figure 2: a and b: One-loop graphs with the insertion of the pion operator c: Tree graph 

with the insertion of the pion operator On,-K{(^)- The operator insertions are denoted by black blobs. 

Let us finally return to the contributions to the nucleon matrix elements from the pure pion 
operators 0^^{a) and 0^^{a). Their chiral dimension is [13] 

D^ = 2L-l + d^ + Y, (dim K(j) " 2) + J]] (d™ Kiv(i) - l) (37) 
j=i i=i 

where > n is the chiral dimension of the pion operator. Because of parity invariance the vector 
operators On.-w couple to 2 or more pions, whereas the axial vector operators 0^^{a) couple to 1 but 
not 2 pions. Starting at order 0{q^~^) the form factors E^^k+i ^md M„ ^ thus receive corrections 
from the one-loop graphs shown in Fig. [2^ and b. For isotriplet pion operators n is odd due to charge 
conjugation invariance. Together with the time reversal invariance constraints on the nucleon form 
factors, one thus finds that the corrections to E^^_^ start at 0{q^) and those to M^^^j^ at order 
0{q), whereas for all other form factors -E^^+i and they are at least of order 0{q^). 

The axial vector operator 0^^{a) contributes to nucleon matrix elements starting with the tree 
level graph in Fig[2j;. The n vectors in the operator are all contracted with derivatives acting on 
the pion field and hence with after evaluation of the graph. The same is true for the corresponding 
one-loop graphs. One thus obtains only contributions to the form factor A£„^„_i, starting at order 
0{q^'^) for the tree graph. With two loops one has graphs where three pions couple to the operator 
on one side and to the nucleon line on the other. Such graphs can contribute to other form factors, 
but only starting at order 0{q'^). 

5 Results for chiral-even isotriplet form factors 

5.1 Relevant operators and graphs 

With the method outlined in Section [4.11 one finds vector operators 
= (^«^)' i^S) rlAT, + . . . , 

<fc+i,o = Ki+i {^'^df^' r^N. - \Mif {iadf td,K [(a5), S^]t^^N, + . . . , 
Oik+1,1 = \Mif {iadf^^ id,N, S'^rtN. + ..., (38) 

where the . . . stand for operators which have fewer total derivatives and as in Section [2] the isospin 
index / = belongs to ^ = and / = 1 to ^ = 1, 2, 3. The axial vector operators are simply obtained 
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by interchanging and t^_, 

= \Mif {iadf^' id,N. S^T^^K + . . . . (39) 



Using the rules (|34l) and the decomposition ([8]) the coefficients in (pSl) and (l39]) are easily identified 
as the tree-level contributions to the form factors E^^k+i^^)^ ^nk ^"^^ ^nk order 0{q^). 

Time reversal invariance implies that E^^^^^-^ is only nonzero for odd k, whereas the other coefficients 
M^^^\ E^^^\ Af^^^"* are only nonzero for even k. For A = we recover the isosinglet operators 



constructed in [T3|. 

According to (j29p insertions of an isotriplet operator with require at least one pion line in the 
graph and hence do not contribute to nucleon matrix elements at tree level. They appear however in 
the one-loop graph shown in Fig. [1}:. When the pion-nucleon vertex in this graph is taken at LO one 
finds zero, because the loop integral is of the form 

d^-^H f40l 

{Iv + w + iO) {P - + iO) ' ^ ' 

whose numerator is proportional to (Sv) = after the integration. Calculating the same graph to 
the next order, one finds that the contribution from the NLO pion-nucleon vertex cancels the one 
with the LO pion-nucleon vertex and an NLO nucleon propagator correction. This holds true for 
all operators with in (I38p and ()39p and only requires that the operators does not introduce any 
dependence on the loop momentum via or u^. 

The operators with r°^_ contribute at tree level and via the loop graphs in Fig. [1^ and b. They 
are constructed such that after the replacement 5^ iA^ they match the structure of the terms in 
the form factor decomposition ([8]). That structure can be changed in loop graphs only when the spin 
vectors in the operator insertion are multiplied by spin vectors from pion-nucleon vertices. This is 
not the case for the graph in Fig. [lb, which originates from the two-pion term in the expansion (1291) 
of Tg_|_ and thus reproduces the spin structure of the operator. Let us show that it is not the case 
either for the graph in Fig. [1^ with LO pion-nucleon vertices. The numerator of the corresponding 
loop integral has the form {Sl)0{Sl), where O contains zero, one or two vectors 5"^ and represents 
the spin structure of the operator. The loop integration turns a tensor . . . Ifj,- into a combination 
of and gxfj_ and thus {Sl)0{Sl) into S^OSp. This preserves the spin structure of O because 

SPSp = ^il-d), SPS'^Sp = ^{d-3)S'^ , SP[S'^,S^]Sp = ^{5-d)[S'^,S^] (41) 

in d dimensions. We also need graphs with one LO and one NLO pion-nucleon vertex, or with two 
LO vertices and an NLO nucleon propagator correction. Restricting ourselves to the terms producing 
the required factors of as discussed in Section we obtain numerators {lv){SA)0{Sl) or 
{IA){SI) O {SI), which give zero after loop integration. 

In summary, the insertion of an operator from (|38p or (|39p into the graphs discussed so far either 
gives zero or contributes only to the same form factor for which it already provides the leading-order 
tree- level result. This is however not true for the graph in Fig [1^ with one NNLO or two NLO 
pion-nucleon interactions. In this case one obtains terms {SA)0{SA) after loop integration, which 
do change the spin structure of O. 
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5.2 Vector form factors 



From Table [D it follows that Ej^ ,^_^^ and M^ f^ can receive contributions from one-loop graphs with 
insertion of operators On,i+i,i-i with / > k and i > 0. With the additional condition (j36p required 
to produce enough factors of A^, we find that the form factors receive corrections of order 0{q^) 
from graphs with LO pion-nucleon vertices and insertion of the operator On,fc+i,0! which already 
gives the tree-level contributions at order 0{q^). By power counting one could also have order 0{q'^) 
contributions from graphs with insertion of On,k+2,-i or On,k+i,-i and pion-nucleon interactions 
at LO or NLO, respectively, but these vanish because the relevant operators come with T^_. The 
one-loop corrections from pion-nucleon operators to E^f}{t) and M^^{t) are then found to be 



M,:?'°' 1 1 



m 



2 

m 
{A-kFY 



m 



(35l + l)log^ + 2gl 



m 



(2^1 + 1) log ^ + 2^1 



+ 0{q^), 
+ 0{q^), 



(42) 



respectively, where we have replaced the bare axial coupling go by its physical value qa as is permissible 
within the precision of our result. Likewise, we could replace the bare pion decay constant F and 
bare pion mass m by their physical values -Ftt and rriT^ (we refrain from doing so for ease of notation). 
The contributions with g'j^ in (j42p are due to the graph in Fig. [1^ and the nucleon wave function 
renormalization, and the contributions without g\ come from the tadpole graph in Fig. [TJa. As in 
|13j we use the renormalization scheme of [22], subtracting 1/e -|- log(47r) -|- ip{2) for each 1/e pole in 
4 — 2e dimensions. 

The form factors -E'^^^^(t) and also receive chiral corrections from loop graphs with pion 

operator insertions. In the notation of [12] the isotriplet operators with lowest chiral dimension ar^ 



n— 3 



fc=o 

even 



where 



LlT'' = U^d^U, 



(44) 



To extract the terms coupling to two pions we use the expansion = id^'7r°_/F + ie'^^'^TT^d^ir'^/F^ + 
0(7r^) and its analog for R'^, obtained by changing the sign of the pion fieldp and obtain 



n— 3 



-2ie'^''^ <^ Ln-i 



{iadT-^Tr^iadir") - 2 ^ K^k {iadf (iaSvr*) {2iad ^-^-'^ {iadir") 



k=0 
even 



ctbcb , 



Using the relations 

4(ia97r^) {iadn^) = {iadfirK^ - tt^ {2iadf vr^ , 

^The normalization of the twist- two operators ((2} used here agrees with the one in [13] and differs from that in [12] 
by a factor of 2. The coefficients bn^k have the same normahzation here and in [12] , 
^We note that the sign of the term with e"*"^ in eq. (32) of [12] is incorrect. 



(45) 



(46) 
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we can rewrite this as 



n-l 



(47) 



fc=o 

even 



where A^^^^ = 2" ^{hn,k — bn,k-2) with 6n,-2 = 0. The coefficients A^^^^ represent the chiral hmit of 
the form factors A'^ ^ (t) which parameterize the moments of the pion GPD as [12] 



/■I " 

dx x^-^Hi{x, t) = Y.m^'Ai^^it) , 



(48) 



fe=0 
even 



where in terms of quark flavors in a vr"*" one has H^.^^ = + H!^ and H^^^ = — H^. Because of 
isospin and charge conjugation symmetry one has 1 = 1 for odd n and / = for even n and therefore 
can omit the isospin index / in ^. 

As discussed after (j37|) . the graphs in Fig. [2^ and b with insertion of 0'^^{a) give rise to corrections 
which start at order 0{q^) for EI^^_^ and at order 0{q) for M^^^^^. Together with (jl2]) and with 
terms due to tree level operator insertions, the complete results to order 0{q'^) read 



1=1(0) 



m 



(47rF)2 



m 



(35l + l)log^ + 25l 



.2 



7=1 (2,m) 2 



7=1 (2,t) 
n,k 



t + 0(g^ 



7=1(0) 



+ ^k,n-l 

where the contributions 



1 



m 



(47rF)2 



m 



(2gl + l)log^ + 25l 



7=1 (2,m) 2 



n,k 



7=1 (2,t) 



(49) 



p7=l(2,7r). X 



odd 



+ 



dr] Tj- 



(2m2 - t) log 



+ 1 



m?'{ri) 



q\ - \)rr?(r\) log 



j=i J-^ 

odd 



m{r]) 



with 



m^(r/) = — — (1 — 77^) 



(50) 



(51) 



are due to graphs with pion operator insertions and LO pion-nucleon vertices. The order 0{q^) 
correction 



«,;-„■<?••'«) 



odd 



+ 



g\ [2w? - t) ( log + 1 ) + (g2 _ 1 _ 4^^^) ^2(^) 



m?{rj) 



(52) 
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is due to graphs with one NLO pion-nucleon vertex or nucleon propagator correction, as well as graphs 
with LO vertices and the subleading part wv^ of the residual nucleon momenta, cf. the discussion 
after ()33p . The terms proportional to g\ in ()50p and ()52p are due to the graph in Fig. [2^, and the 
other terms to the graph in Fig. [2)3. Our expressions (|49p and (j50p agree with the results in , 
where the order 0{q^) corrections to E^^_^ and the order 0{q) corrections to M^^^^^ are given. 



5.3 Axial form factors 

Using Table [Hand the condition (j36p . one readily finds that the chiral corrections of order 0{q'^) to the 
form factor ^ are obtained from graphs with LO vertices and insertion of the operator On,k^ 



1,-1) 



which already gives the tree- level contributions at order 0{q^). Together with higher-order tree level 
insertions we get 



/=1(0) 
n,k 



1 



m 



(47rF)2 



m 



{2g\ + 1) log ^ + g\ 



+ E., 



7=1 (2,m) 2 



n,k 



+ E, 



I=l(2,t) 
n,k 



t + 0{q^ 



(53) 



in agreement with [14]. The discussion of contributions to is more involved; for the isotriplet 

case it proceeds in close analogy to the isosinglet case analyzed in [13] . According to Table [1] and the 
condition (j36p . one obtains order 0{(f') corrections from graphs with insertion of On,k+i,i and LO 
vertices. Further corrections are due to graphs with insertion of On,k+i,-i and two NLO pion-nucleon 
vertices or nucleon propagator corrections, or with one NNLO pion-nucleon vertex generated by (j35p . 
Graphs with the same operator insertion and two loops or one loop and a pion propagator correction 
from could contribute by power counting but do not produce the required factors of (see 



[13]). Graphs with insertion of On^k+2,0 or On,fc+i,o and pion-nucleon interactions at LO or NLO give 
zero because these operators involve r^_, as discussed after (|40]l . Graphs involving On,fc+2,-i and 
NLO pion-nucleon interactions do not contribute to ^ due to time reversal invariance, since the 
operator is only nonzero for odd k and the form factor only for even k. Finally, graphs with insertion 
of (Dn,k+3,-i and LO pion-nucleon vertices contribute to E^ but not to ^ as discussed at the 
end of Section 15. li 

Together with higher-order tree-level insertions, the one-loop graphs with On,k+i,i or On,k+i,-i 
thus give the full result at order 0{q'^) for the form factors with A; < n — 1, 



m 



(47rF)2 



m 



{2g\ + 1) log —+g\ 
1^ 



+ E, 



%7 log — ?r + M , ^ ' m + M , t + Oiq 

3(47ri<')2 /i^ ' ' 



n,k 



(54) 



The form factors M^^_^ require a separate discussion because they receive a contribution starting 
at order 0{q~'^) from the one-pion exchange graph in Fig. [2]:, as discussed at the end of Section W% 
The relevant operator is given by 



1 + 



("2^2 



+ c„ Tr x^ 



+ 0{q 



(55) 



with odd n and 



(56) 
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where 6n,n-i is the same as in (I43p because of parity invariance. is the renormahzed low-energy 
constant from the pion Lagrangian ()17p and appears in the expression of the axial current, 



^QitMibT q = 



1 + Jl^XH 



+ Oiq'), (57) 



so that 6i^o = 1 ^-nd ci = 0. One can readily derive ()57p by coupling the Lagrangian to an external 
isovector axial field as usual [12], which implies = i{v)d^u — udu^^ + u^a^u + ua^u^ . As an 
aside we note that the correction with in the axial current (j57p would be different if one used the 
pion Lagrangian from [22|, where the term involving this low-energy constant reads 

^Tr'- - — ^ 

4 

In the present work we follow Ref. [18] and use the Lagrangian (jl7p from [T7]. It differs from (j58p 
by total derivative terms and terms that vanish by the equation of motion. With the full Lagrangian 
given by C-,^ + C-j^i^ the equation of motion for the pion field involves terms bilinear in the nucleon 
fieldS so that a change of C-,^ using the equation of motion induces a corresponding change in Ct^n- 
We also refer to the discussion in [23] . 

The coefficients hn,n-i and c„ in (j55p appear in the moments of the twist-two pion distribution 
amplitude (ji-wix), which are defined by 



/ 



and 



so that 



^ e*-'?M(^f (p)| q^-r^a) ^j^r^qiva) |0> = -iS'^'^Mx) (59) 



K = 2"" / dxx"-^^^{x), (60) 



{Ap)\O^^Ja)\0) = -2i5^\aprF^Bl . (61) 



Calculating the leading chiral correction to this matrix element one finds that l\ only appears in the 
expression of the pion decay constant, 

F^ = f{i- log ^ + ^ ^I(^) + ' (62) 

and Cn only in the correction to the moments 

Sn= Vn-i(l + 4"i'c„)+0(m^). (63) 

The definition of the pion decay constant implies = 1 to all orders in the chiral expansion. 

Returning to the nucleon form factors M^^^_^, one readily finds their lowest-order contribution to 

M-<-'(*) = B„"'»'i^ (64) 

with Bn^^^ = bn,n-i- The leading corrections to this come from a number of higher-order operator 
insertions and loop graphs. A tadpole insertion in the pion line of the graph in Fig. [2]: and the 



*With an arbitrary matrix X in isospin space and X — X — | TrX, the leading-order equation of motion reads 
2i Tr [V^, ^t"] X + Trx-X = F'^ iV„ [{vu) + 4igo (SV), X]N^ ^ 4i(?o-F"^ d" (7V„ S^XN^) . 
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propagator correction from the ^3 term in the pion Lagrangian ()17p result in a shift 
where 



ml in 



m 



1 + 



2(47rF)2 



(65) 



The pion propagator correction from the I4 term in (jl7|) cancels against the ^4 term from the operator 
(I55p . A tadpole insertion at the operator vertex gives a chiral logarithm as in (1621) . Further corrections 

(3) 

are due to loop corrections to the pion-nucleon vertex, to tree- level insertions from C^^, and to the 
factors A/"^ and Zj\f in the matching formula ()2ip . Together with the tree-level insertion of the pion- 
nucleon operator „ from (j39|) we obtain 



4M250 



mi 



(47rF)2 



m 



(25o' + l)log^+5o' 



+ 2^250-^ (2rf 



16 ~ '^is) 



Sm'd^s + 4m^c„ ) + M^'^.^ + 0(g) (66) 



7=1(0) 



with the low-energy constants (i^g(/i), d^^ and d2^{^) from [21j. To make the pion mass dependence 
fully explicit, one should replace = Mq — Sm^MoCi + O(m^). Conversely, we can use (j63p and 
the one-loop expression 



m 



(47rF)2 



(25^ + 1) log ^ + 5o 



+ Arr?gQ^dr^f^ - Sm^dT^^ + 0{m^) (67) 



of the axial coupling (see e.g. [24J) to rewrite the result as 



Mi=Li(t) = Bl (1 - 2mlgfd,,) + mItS^ + 0{q) 



mt, 



t 



(68) 



r/=l(0) 



in terms of the physical quantities B^, qa, M, m,r and the low-energy constants dig and M^~^_^ . 
With the transformation (jlOp and the definitions ([5]) and (j60p of moments one finds 

'""^ ^^{l-2mlg^%,)+E^='('\x,0 + 0{g), (69) 



K\x\<0 



2C --y^J mi-t 

which generalizes the well-known relation from j25tl26j to next-to- leading order in the chiral expansion. 
For n = 1 our result (|68p is consistent with the one for the pseudoscalar form factor in |18[ [27] . We 
also agree with the result of Ref. if in their eq. (66) one adds a term 5m,2k | gAM'^{r'j^ {z^^) j 2"^^ . 



6 Chiral-odd generalized parton distributions 



In this section we consider the general parton distributions associated with chiral-odd operators. As in 
the previous sections we restrict ourselves to the twist-two sector. The relevant GPDs of the nucleon 
are defined bjH 



^ e"''("^)(Ar,(p')| g(-^r?a) ftAa^i^"'^ A(M I^J-(^^)> 



2aP 



bxa^ u{p') 



2M 



2M2 



M 



u{p) . (70) 



^We have traded the distribution Et in the original decomposition _28_ for the combination Et = Et + 2Ht, which 
naturally appears when representing the distributions at ^ = in terms of densities in the impact parameter plane [8]. 
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In addition to the light-like vector a we have introduced a vector b satisfying ab = 0, and for brevity 
we have suppressed the arguments x, t of the distributions]! Their Mellin moments are related to 
the matrix elements of the chiral-odd twist-two operators 



A/ii.../t„ A/^i /ii.../t„ 



(71) 



where T and S are defined as in Section [2] and where A denotes antisymmetrization, A;^^ t^^ = 
i^^A/i_^/iA^_ These operations are conveniently implemented by contraction with the auxiliary vectors 



a and 6, given that for any tensor t^i^i---t^" satisfying t^i^^^^' 
bxa^,a^2.--a^„ T A S t^'^i-^" 

A/Xl.../^n A/Xl /il.../in 



^ n n 

'^^ i=2 i=2 



n + l 
2n 



(72) 



where symmetrization in /i2 ■ • • A*n is guaranteed by contraction with identical vectors, and where trace 
subtraction terms are removed by the conditions a? = ab = 0. The (n — 1) terms of the second sum 
give zero due to the antisymmetry of t in its first two indices. We therefore define the contracted 
operator 



In 



Jb^a''iaxf,{iaDy 



(73) 



whose nucleon matrix elements are parameterized by 

n-l 

{N,{p')\ OUa,b) \N,{p)) = ^^(aP)"-^-! {a^f bxa,u{p') 

k=0 



r 



2M 



B 



Tn,k 



ia^"" A^A^ - iaf"" A^A^ ~j 

2M2 ^"•'^ + 

The moments of the chiral-odd GPDs are then expressed as [29] 



^A ptj. _ pX^ii 

B' 



M 



Tn,k 



u{p) . (74) 



dx X"-1 Ht{x, e, t) = ^(-20^^ ATn,k{t) , / dx Et{x, t) = J2(-'^0^ BTn,k{t) , 

fc=0 -^"^ fe=0 
even even 

/I n— 1 n— 1 

(ixx'^-^^Hx.e,*) = E(-20'^Tn,fc(t), / dx:E"-i^T(x,e,t) = V(-20'5Tn.fe(t), (75) 

-1 tTn J-1 t~, 



k=0 
even 



k=l 
odd 



where we have omitted isospin indices I in the distributions and form factors for ease of writing. The 
restrictions to even or odd k for the form factors reflect that Ht, Et and Ht are even in ^ and Et 
is odd in ^ due of time reversal invariance 
Using the relations (9) in [13] and 



u{p)ia^'^u{p) = Uv{p) 



27[S\5^] + 



2M 



+ 



2M2(1 + 7) 



Uvip) (76) 



^Instead of contracting a^'^ with auxiliary vectors one often takes definite indices a'"'', where i = 1,2 denotes a 
transverse component and + the plus-component in light-cone coordinates, i.e. a*"'' = (ct'" + (j'^)/\/2. 
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one can rewrite the decomposition ()74p in terms of the heavy-baryon spinors ([7]) and obtains 



n-l 



{Niip')\ OUa,h) \N,{p)) = Y.^M^)'''''' ^''^)^'''~' («^)' ^^a^. 



fc=0 



[5^(SA)]A'^-[5^(5A)]A^^, 



2M2 



Tn,fc 



[g^(5A)]^;/--[g^(gA)]^;^ , 7; 

■ 7 -^Tn.fc + 



2M 



■Mr 



Tn.k 



Mp) (77) 



with new form factors given by 



Tn,k 



B. 



Tn,k ' 



MTn,k = (1 + 7) ^Tn,k + BTn,k " '^^Tn,k 



(78) 



or equivalently 



B 

A 



Tn,k 



Tn,k 



1 
1 

2^ 



A2 ~ ~ ■ 



^Tn,k ~ ^Tn,k 



TT^^Tn,k 



(79) 



We finish this section by defining chiral-odd GPDs in the pion, 



dr] 
47r 



,ixri{aP) / c( / 



2aP 



(80) 



where as in the nucleon case, isospin 7 = corresponds to A = and / = 1 to ^ = 1, 2, 3. In terms 
of quark flavors in a tt~^ one has = E^^ + E^^ and E^^ = E^^ — E^^, with the definition 



2aP ' m. 
and its analog for d quarks. For the local twist-two operators one has 



(81) 



{7T%p')\ OUa,b) \7r\p)) = iTr(r^rV) b,a. 



with 



n-l 



Y,{aPT-'-\a^f BJ^^^,{t) (82) 



fc=0 
even 



/ dxx--^EUx,i,t) = X;(20'Sfn,.(t), 



(83) 



k=0 
even 



where the restriction to even k is a consequence of time reversal symmetry. Due to isospin and charge 
conjugation invariance, n is even for I = and odd for / = 1, so that we do not need an isospin label 
for B^l. 
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7 Chiral-odd effective operators 

In this section we explain how to construct the operators in the effective theory that match the chiral- 
odd quark operators (ITT]) , closely following the strategy used in Section HI To this end we first match 
the operators 



with open isospin indices i, j, which involve quarks of definite chirality. The corresponding uncon- 
tracted operators \qjia\^j_^{\ ± 75) iD^^ . . . iD^^ qi do not have definite twist, but according to ([72|) 
their twist-two part is readily projected out in bxa^O^j^ ,^{a) and h\a^O^^^{a). The operators ([8^ 
transform as 

- VRO'n'W^Wl 0^'^,Ja) - V,Ol^^^^{a)vl (85) 

unter chiral rotations and are transformed into each other by parity. Because o"^^75 = — ^ie^^a/j o""'^ 



they obey the duality relations 

<,„(«) = <,n(«) oI'rM = Of^Ja) (86) 

The operators 0^^{a, h) from ([75]) . which correspond to twist- two and to definite isospin, are obtained 
as 

O^Ja, 6) = bxa, Qi'^^{a) , Qi'^^^{a) = TVr^ (O^'^^Ja) + Ol%^{a)] . (87) 

They will involve the combinations 

= uV^ J ± ut\ , (88) 

whose expansion in pion fields reads 

= 2 - TT'^vr'^ + ©(vr^) , = -| vrV^ + 0(^3) ^ 

= 2r'^ - vrW + ©(vr^) , = -| + ©(vr^) . (89) 

As for the chiral-even case discussed in Section 14.11 the operators which match (j84p in the effective 
theory and contribute to nucleon matrix elements are either bilinear in the nucleon field or contain 
only pion operators. We treat the two cases in the following two subsections. 

7.1 Pion- nucleon operators 

The effective operators which are bilinear in the nucleon field and transform as (|85p can be written 
in the form 

(oXnia))., = (N.O.u)^ (uO^N,)^ , («)),, = (NvOW), {n^O'.N,)^ , (90) 

where Oi, O2 involve the fields u^, x± covariant derivatives and transform like under chiral 
rotations. 0[ and O2 are related to Oi and O2 by parity. One can rearrange the covariant derivatives 
in O^j^ ^ and O^'^ ^ such that they act either as total derivatives 9^ or in the antisymmetric form 

v^ = i(v^-v^). 
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Table 2: Overview of contributions to the chiral-odd form factors. As in Tabled] the restriction in the 
second column is due to time reversal invariance. A^a is the number of factors (aA), (bA) and (SA) 
in the decomposition ([77|) . One must have /> A; — 1 for ETn,k, MTn,k and / > for ETn,k, MTn,k, 
and i > for all cases. The corresponding graphs contribute to the form factor at order 0{q'^) with 
d > Di^i - Na and A,j from (f33]l . 



form factor k N^. operators 

ETn,k 
MTn,k 
ETn,k 
MTn,k 



Odd k + l Q'^l^ 

even k + 1 . 

even k Q^^i^ 

even k + 2 Qt^,. 



To obtain the general form of ^ and „ it is sufficient to construct corresponding operators 
O^^ that involve no e tensor and either no spin vector or two spin vectors in the form [S\ , 5^] . 
Operators with one e tensor and one spin vector can be brought into this form by using the third 
relation in (j30p and rewriting the resulting product of two e tensors in terms of products of metric 
tensors. Terms in O^'^ „ and with an odd total number of e tensors and spin vectors are 

then readily obtained by adding the dual operators ^ie^^^ap^"^ with coefficients determined by the 
relations ([86l) . using tlicit "^^C^^Q-^ ^ic^^ any antisymmetric tensor t^^ . 

Following the procedure of Section 14.11 we decompose the pion-nucleon part of the operators 
Qn'^'^ia) as 

n-l 
k=0 

where we have omitted superscripts A for ease of writing. The operator Q^'^(a) is the contraction 
of a tensor of rank k + 2 with k vectors a and may not contain any factors (av). The minimal 
number of vectors dp, Vp, Up in Q^^j^{a) is k — 1 and must be accompanied either by the tensor 
v^lS^^, (aS)] — vf^[S^, {aS)] or by its dual ie-^t^'^l^ VaiSp, (aS)]. In the first case one obtains however the 
structure {av)[{bS), (aS)] after contraction with bxQf^, which also appears in bxa^ {av)""'^ Q^'^k-ii^)- 
An analogous statement holds of course in the case of the dual tensor. We can therefore restrict 
ourselves to operators Q^^i.{o) with at least k vectors dp, Vp, Up, and thus further decompose 

oo 

Qtk{^) = Y.^i-'QtkM)^ (92) 

i=0 

where Q^^iio) has chiral dimension k + i. The power counting for graphs with a certain operator 
insertion proceeds in close analogy to Section 14.21 and is summarized in Table O Comparing the 
number of factors {av) in (j9ip and in the decomposition (j77p . one obtains the restriction I > k for 
the operators Q^^ j(a) that can contribute to E^n^k and Mxn,k- For Exn,k and M^n^k the restriction 
is / > A; — 1, where the case I = k — 1 requires that the graphs with insertion of Q^nicL) produce no 
factors of or v^. 
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7.2 Pure pion operators 



Pionic operators which transform according to (185^ can be written as 

0Xni<^)=n0u, Ol%Ja) = u^O'u^, (93) 

where O and O' are related by a parity transformation, transform hke under chiral rotations, and 
are constructed from the fields Up. x± and covariant derivatives. We can restrict the derivatives to 
act only on fields inside O and With the duality relations ()86p one finds that the pure pion part 
of the operator Qn''^^{cL) can be brought into the form 



Tr 



T^+V^^'ia)] + e^^°^ Tr fr„lK./3(a)j or TV [t^ A^t'ia)] + e^^°^ Tr \t^^ A„^(a)j , (94) 



2 



where Vx/i and A\fj_ respectively behave as a tensor or a pseudotensor under parity and are constructed 
from Up, Vp and x± without the e tensor. One readily finds that the terms without e in (j94p couple 
to an even number and the terms with e to an odd number of pion fields. Vxp{a) and A\pi^{a) are 
tensors of rank n + 1 contracted with n — 1 vectors a. In the following we consider the terms with 
the lowest chiral dimension in the pure pion part of 0^^{a, b). These terms contain no fields x± and 
have the vector indices of all n + 1 factors Up or Vp contracted with either a or b. 

To calculate matrix elements of these operators between two nucleons or two pions at one-loop 
accuracy, we only need terms that couple to at most four pions. Terms coupling to three or four pions 
can appear in tadpole graphs. Such graphs are only nonzero if the pion fields in the operator which 
couple to the loop have no derivatives acting on them. This is because the corresponding loop integral 
has a numerator of the form Ip^ . . . Ip^ , where I is the loop momentum. After the loop integration, 
one obtains zero for odd m and for even m one obtains a combination of metric tensors, which gives 
zero when the vector indices are contracted with a or b. 

Since the derivatives with indices A and /i in the antisymmetric tensor Qn'^^{a) cannot act on 
the same pion field, one readily finds that operators coupling to one or three pions do not contribute 
to matrix elements between two nucleons or two pions. For the same reason such operators decouple 
from matrix elements between the vacuum and a single pion, which reflects the fact that there are no 
chiral-odd pion distribution amplitudes of twist two. 

It remains to construct operators Vxp{a) and Axp{a) from Up and Vp, which must have at least 
one factor Up because the covariant derivatives must act on some field to give nonzero, and less than 
three such factors because of the restriction just discussed. The operators with one factor of Up are 
of the fornjl 

{aVf^V^{aVf^u^'-{\^^l) or {aVf^V^{aVf^V^'{aVf^au)-{\^^JL). (95) 

In both cases we can use the commutator identity [V", V^] O = \{[u'^ ,u^]0 - O [u", li^]) to bring 
the vectors with indices A and next to each other. The commutator terms do not contribute to the 
matrix element in question since they involve three or more vectors Up. The remaining term involves 
either V^u^ — V^u^ = or [V'^, V^] . . . (au) and thus do not contribute either. 

The only relevant operators contain hence two vectors Up. According to our above discussion, the 
Fp part of any factor Vp does not contribute in this case, and the derivative must act on the pion 
fields in Up which already carry a derivative. For the matrix elements in question, V^ (au) is hence 
equivalent to (aV) u^. The same holds of course for the index We thus find that the operators of 



^Other terms can be brought into this form using identities such as dp{uOu) — u([Vp,C'] — ^{up, O}) u. 
®For simphcity we write from now on V^O instead of [V^, O] if O transforms hke under chiral rotations. 
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lowest chiral dimension can be written as 

Qn:''(^) = ^ E ^Tn,k Tr (^aV)'= i;^^ + . . . (96) 

A:=0 
even 

with 

Y^"^ = (2mV)"-^~i - u^" (2iaV)"-'=-i , (97) 

where the . . . denote terms not contributing to two-nucleon or two-pion matrix elements at tree level 
or one loop. We note that the coefficients hTn,k have nonzero mass dimension and are of order (47ri<')~^ 
in the sense of chiral power counting. They give the tree-level contribution at order 0{(p) to the pion 
form factors ^ (t) defined in (f82l) , 

^?S = (-ir+'2"-'="'"^-W' (98) 

where n is even in the isosinglet and odd in the isotriplet case. The restriction to even k in (j96p 
corresponds to the one in ([52]) . 

We can now apply the power counting formula (j37p with ciTr = n + 1 to the operators just con- 
structed. Taking into account the restrictions of even or odd n or k for the different form factors, we 
find that the corrections from pion operator insertions start at order 0{q) for M^^^_-^ and at order 

0{q^) for M^^\_^, ^Tn^n-i ^"^^ -^^'^ Other form factors they start at order 0{q^) or 

higher. 

8 Results for chiral-odd form factors 

Using the construction described in Section [71 we find pion-nucleon operators 
Qtifi = ETnl^iadf (N,[S\S^^]T^^N, + K[v^S^-v>'S^)TtN.)+... , 

+ (^a^)' dp(v>'N,[S\ SP]t^_,N, - v^N,[S>', ^^Jr^+iV, + {terms with r^}) + . . . , 

<fc,2 = -\M^nl dp(^d^'N,[S\ S^lr^+iV, - 9^iV45^ 5^]r^+iV, + {terms with r^}) + . . . , 

(99) 

where the . . . denote terms with a smaller number of total derivatives. The coefficients in (I99p are the 
tree- level contributions at order 0{q^) to the respective form factors and therefore only nonzero for 
even or odd k as given in Table [2j The terms with in the last two lines of (I99|) are rather lengthy 
and not given here. Indeed, one finds that none of the operators with in (j99p contributes in one- 
loop graphs with pion-nucleon interactions at LO or at NLO. Such graphs have the form of Fig. [It and 
give zero for the same reasons discussed after (|30|) for the case of operators with r^_. The discussion 
at the end of Section [5.11 also applies to the operators with r^^ in ()99p . so that their insertion into 
graphs with LO pion-nucleon vertices or with one NLO pion-nucleon interaction contributes only to 
those form factors for which they already provide the tree- level result at order 0{q^). 

An operator Q^'^ ■ in ([99|) has at most / -|- i partial derivatives, so that the condition ([36|) holds 
also in the chiral-odd case. Together with the power counting following from Table [21 one again finds 
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that one-loop corrections to all form factors start at order 0{q^). One finds that the order 0{q^) 
corrections to MTn,k and ETn,k come from Qn,k,i, whereas those to ETn,k and MTn,k come from Qn,k,o 
and Qn,k,2j respectively, with pion-nucleon interactions taken at LO in all cases. Additional order 
0{q^) contributions to M^n^k come from graphs with Qn,k,o and two pion-nucleon interactions at NLO 
or one pion-nucleon interaction at NNLO (only the t^_^_ part of the operator is found to contribute). 
Contributions from the or MTn,k+i are possible by power counting but turn 

out to be zero. Other contributions at order 0{q^) which are possible by power counting involve 
at most one pion-nucleon interaction at NLO and do not appear for the reason given at the end of 
the preceding paragraph: there is no correction to Mxn.k or Exn,k from Qn,k+i,Oj Qn,k,o, Qn,k~i,i or 
Qn,k-i,2 and no correction to MTn,k from Qn,k+i,o, Qn,k+2,o, Qn,k,i or Qn,k+i,i- 

Taking into account the graphs with pion operator insertions shown in Fig. [2^ and b, we finally 
find 
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and 
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where for brevity we have written ... to denote analytic terms proportional to m'^ or t and corrections 
of order 0{q^). The analytic terms are due to higher-order tree-level insertions as specified below 
(I34p . The contributions from pion operator insertions read 
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with n even in the isosinglet case and 
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with n odd in the isotriplet sector. As remarked in [14j, the corrections from pion operator insertions 
are very similar for the chiral-odd and chiral-even form factors. We find a correspondence 
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for the terms in (jl02p and (jl03p when interchanging Mbj,^ <-» 2 ^jA^J\^_y 

Let us compare our results (jlOOp to (jlOSp to those in Ref. [m] , which gives the corrections of order 
0{q) for Af/^^^ and of order 0{q^) for all other form factorslfl We agree with the expressions given 
there, except for the corrections from nucleon operator insertions without a factor g\ in the isosinglet 
form factors, which are absent in [Ijj, and for the corresponding term in M^^^^, where we have a 
different coefficient. These corrections are due to the tadpole graph in Fig. [T]3, with the pion-nucleon 
vertex generated by the two-pion terms in the expansion (|89p of t^j^ and r"_^. Since this vertex has 
no spin or momentum structure, the corresponding corrections must be the same for all form factors 
with a given isospin. 

Let us finally give the corrections of order 0{q^) to the chiral-odd GPDs of the pion. They are 
given by the one-loop graphs shown in Figure [3] with insertion of the pion operators (j96p and from 
tree-level insertions of operators with chiral dimension n + 3. For the form factors (I82p we find 
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for even n, and 
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®The tensor form factors in (14] are related to those introduced here by Ai'^ ~ Et, ~ Mt jl, 
up to terms suppressed by factors of order /M^) and by = Et- 
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Figure 3: One- loop graphs contributing to two-pion matrix elements of the pion operator Qn'^^{a) 
in (j96p . The operator insertion is denoted by a black blob. 



for odd n, where the . . . stand for analytic terms from tree- level graphs and for corrections of order 
0{q^). The corrections going with log / iJ? are due to the tadpole graph in Fig. [3^ and are inde- 
pendent of k. The term involving an integral over r] is due to the graph of Fig. [8)3 and can only occur 
for k = n — 1 (and thus only in the isotriplet case). This is because the operator insertion on the pion 
line cannot produce any factor (aP) and the four-pion vertex can only produce one such factor after 
the loop integration. 



9 Results for moments of nucleon and pion GPDs 

In this section, we rewrite our results in terms of the form factors An^k, Bn^k, Cn, ^n,fc) Bn^k and ATn,k, 
BTn,k: ATn,k: Bxn,ki which describe the moments of GPDs in commonly used parameterizations. We 
give expressions for the value and the first derivative of each form factor at t = 0, which should 
be useful for applications in lattice QCD. The corrections obtained from graphs with pion-nucleon 
operator insertions are completely specified in this way, because they are independent of t. In the 
following we will use the abbreviation = 4ttF. 

For a convenient overview of results we also reproduce the expressions for isosinglet distributions 
from [13] here. Together with ([9]) , (|10p and the expressions in Section [5l we find that up to corrections 
of order O(m^) the chiral-even vector form factors at t = have the form 
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The labeling of coefficients with superscripts (0), (2,m), (l,7r) and (2,7r) follows the same pattern as 
in Sections [5] and [8j The contributions from pion operator insertions read 
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where n is even, and 
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where n is odd. The A^^^\j are the leading terms in the chiral expansion of the pion form factors in 
M and fulfill the relations [HI 
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where Bn^^^ is the n-th moment of the pion distribution amplitude to leading order in the chiral 
expansion, as introduced in Section 15.31 Estimates for the values of the low-energy constants q 
appearing in (jlOSp and (|109p can be found in [30] . 
For the axial form factors we have 
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with corrections of order O(m^), and 
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Note that the implicit pion mass dependence from B^, M, g^ and rriT^ is relevant within the accuracy 
of this expression. Numerical estimates of the low-energy constant dig are given in [31]. The derivative 
of B^^_i{t) at t = reads 



d,Bi='_,iO) = Bl (1 - 2mlg-j,'d,,) + 0{m-^] 



(113) 



where the order 0{m ) corrections are due to terms of the form 0{q )/ {m^ — t) in B^^_i{t). Using 
Bf = 1, we obtain a ratio 
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which involves only physical matrix elements and is independent of any low-energy constants. It 
would be interesting to test this relation in lattice QCD calculations, as this would indicate how well 
the chiral expansion works at a given pion mass. 

The derivatives at t = of the remaining chiral-even form factors have nonanalytic contributions 
in the pion mass only for 
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with n even, and 
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with n odd. All other chiral-even nucleon form factors receive only corrections from pion-nucleon 
operators, so that their derivatives at t = are given by the appropriate coefficients with superscript 
(2, t), which are due to tree- level contributions. 

For the chiral-odd nucleon form factors at i = we find 
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with n even, and 
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with n odd. Our results for ^^"^^(0), ^^"^^(0) and Alf^Q{{)) reproduce the expressions in [32] for the 
distributions of unpolarized, longitudinally and transversely polarized quarks and antiquarks in the 
nucleon. The derivatives at t = of the following form factors have nonanalytic contributions in the 
pion mass: 
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where corrections are of order 0{m) and 
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with even n and 
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with odd n. As a consequence of the relations (|lU4p . we find the following correspondence between 
the corrections (|119p . (jl20p . (jl22p . (jl23p from pion loop insertions to chiral-odd form factors and 
their chiral-even counterparts (fTOHD . (fTn9]) . ([TTH]) . (fTTTD : 
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when the low-energy constants are interchanged as Mhj,^^_j <-> jA^^^_-. 

Let us also give the expressions of form factors and their derivatives at t = for the moments of 
pion GPDs. For the chiral-even moments, the expressions given in [12] result in 
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with corrections of order O(m^) and 
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2m? , 



A2 

X 



i=2 

even 
7r{0) 



-.7 J ^4^(0) 



odd 



(126) 
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Using the relation (jllOp one thus has 



<n-l(0)=<n-l + Tylog-3 



, ^7r(2,m) 2 I /n/' 4 



(127) 



with n odd. For the chiral-odd moments we have with (jlOSh and (jl06p 



„^(0) / -I 3m^ 



^-(0) , ^ 



/i2 J ^ Tn,k ^ ^ "k,n-l^Tn,n-l 



2AI 



for even n, 
for odd n, 



(128) 



where corrections are again of order O(m^) and 



rn,n-l 



A2 

X 



odd 



(129) 



The only nonanalytic contributions in the pion mass for the derivatives of form factors are 



X j=2 
even 



diA„,„(U) - + ^2 (^lOg ^2 + i j ^ ^- + 1 ^n,n-, + i2A2 L./ (j + l)(j + 3) ^^'^-^ 

even 

1 / 2 \ 1 

d.Ain-m = a:j^ - (log ^ + 1) E 2--^ , 

odd 

1 / 2 \ " 1 



.7r(0) 
^Tn,n-j ' 



odd 



where the second index is always even and corrections are of order 0{m?). 

Let us finally take a look at moments of parton distributions whose values are fixed by quantum 
numbers for arbitrary values of the pion mass, see e.g. |321 [33] • For yl]^Q(0) one readily finds that 
the explicit chiral logarithm in (jl25p cancels against the one in ()126p . This is required to ensure the 
quark number sum rule 



^lo(O) 



dx (n^ — — + d^) = 2 



(131) 



in the pion, which also implies A^^^^ = 2 and jd^Q^'*"^ 
(|lU7p is consistent with the quark number sum rules 



0. With this one also finds that our result 



A{J'{0)=[ dx{u-u + d-d) =3, A{-^\0) = [ dx {u - u - d + d) = 1 , (132) 
Jo ' Jo 

in the proton, provided that A\~q^^^^ = 3, A\~q^^^^ = 1 and A^'^^^'^'"^^ = A^'^^^'^'"^^ = 0. 

10 Summary 

In this paper and its companion [13] we have calculated the chiral corrections to the full set of twist- 
two generalized parton distributions in the nucleon, using heavy-baryon chiral perturbation theory. 
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For each form factor parameterizing the moments of these distributions, our results include the order 
0{q'^) relative to its lowest-order expression. We have presented a detailed account of the power 
counting and of the operators that can contribute to the chiral order we consider. We find that the 
operator structure is relatively simple in the basis of form factors specified by ([9]) and (|78p . With 
the exception of M„ and Mnpn^k only those pion-nucleon operator insertions contribute to the loop 
corrections of a given form factor which already provide its lowest-order expression at tree-level. 
Furthermore, only operators with or from (j28p and (|88p contribute, but not those with 
or T^_. Our analysis also shows that these simplifications will no longer hold at higher orders in the 
chiral expansion. 

Expressing our results in the basis of form factors parameterizing the moments of the usual nucleon 
GPDs, we find that with the exception of A^^^ and C^^'' all form factors receive chiral corrections 
from loop graphs with nucleon operator insertions (see Fig. [1]). They are of relative order O(g^) and 
contain logarithmic terms m?' log(m^//i^), but are independent of t and of the indices n,k. In several 
cases these corrections involve a mixing between different form factors: i?ri,fc receives corrections 
involving not only its own lowest-order expression but also the one of An^k-, as seen in (I107p . Likewise, 
there are corrections to i?n,fe from An^k) to Bxn,k from Afn^^, and to ^Tn,fc from Axn,k + ^B'pn,k- We 
note that no such mixing occurs for the linear combinations ^n,fc + Bn^k and Axn,k + Bj'n,k- 

Further corrections are due to loop graphs with pion operator insertions (see Fig. [2^ and b). They 
only occur for form factors which are accompanied by the maximal number of vectors in the 
decomposition of the associated matrix element, or by one factor less. Due to the quantum number 
restrictions for pion operators, they only occur for even n in the isosinglet and for odd n in the 
isotriplet sector. Corrections starting at order 0{q) are obtained for C^^^, B^^_i and A^\_^, and 

corrections starting at order O(g^) for B^=J^_2, ^Tn°n-2' -^Tn°n-i' ^n^i-i ^rZ^n-i- order 0{q^), 
the corrections for C^^^ involve the low-energy constants ci, C2, C3 from the pion-nucleon Lagrangian 
(|18p . whereas those for -B^^^x and Aif^^_^ involve C4. The corrections from pion operator insertions 
depend on t. They are responsible for a nonanalytic pion mass dependence of the derivatives of 
form factors at t = 0, namely a 1/m behavior for (?jC^=°(0), df-B^'^_^{0) and ^^A^f^\_-^^{0) and a 
log(m^/^^) behavior in the other cases. We note that these corrections also determine the onset of 
the two-pion cut at timelike t for the form factors in question. 

The pseudoscalar form factors fi^^Li receive corrections from one-pion exchange (see Fig. [2]:). 
They take the very simple form (j68p when expressed in terms of physical quantities. In particular, 
we find that the ratio dfB^^_i{0) /5ji?{o^(0) of derivatives is given by the moment B^ of the pion 
distribution amplitude, with corrections of order m^. It would be interesting to test this prediction 
of chiral symmetry in lattice QCD calculations. 

We have finally evaluated the corrections to the chiral-odd pion GPDs at order O(g^), thus com- 
plementing the calculation for the chiral-even sector. A compilation of our results for the values 
and derivatives at t = of all moments of nucleon and pion GPDs is given in Section [9l 
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